THE INVERSE PROBLEM FOR THE LOCAL RAY TRANSFORM 

Cn ! HANMING ZHOU 

o 

*ivj ' Abstract. In this paper we consider the local X-ray transform for general 

' flows. We extend the results on the local and global invertibility of the geodesic 

i^ ' ray transform proved by Uhlmann and Vasy |13l to the X-ray transform for a 

^3( , general flow. The key improvement is that our argument for the ellipticity of 

'^1>( 1 the conjugated operator Ap (which is dcflned in Section 2) can be applied to 

I -^ ' flows other than the geodesic flow. 

(N 

o 

^"j ' Let X be a domain in a Riemannian manifold {X,g) of dimension > 3. In this 

rS^ ', paper we consider the local inverse problem for the X-ray transform for a general 

C^ ' flow. A general flow (j)t = (7, 7') on SX (the unit sphere bundle on X) satisfies the 
following equation 

(1) ^-^(^'^')' 

K^ ■ where 7 is parameterized by arclength, E{x, v) e T^X is smooth on SX. We call 

-^ , the curve 7 on X that satisfies equation ([T]) an £ -geodesic. Thus we also call the 

f^ ' X-ray transform for a general flow the £ -geodesic transform^ for an ^-geodesic 7, it 

t~^ , is defined as 



^ 



1. Introduction 



^ . If{l) = filit))dt. 

o 

CO ' We extend E onto TX\o (the tangent bundle excluding the zero vector for each 

fiber) by defining E{x,v) = \v\'iE{x, -A-) for {x,v) <E TX\o, we still denote the 
extended one by E. So E is homogeneous of degree 2 in v. When E — 0, this gives 

k> ' the usual geodesic flows. Let x G dX, we say X is £-convex at x if 

A{x,v) > {E{x,v),v{x)) 

for all V e Tx{dX)\o^ where A is the second fundamental form of dM , v{x) is the 
inward unit vector normal to dM at x. If the inequality is strict, then we say X 
is strictly £-convex at x. (These terms are not standard, we only use them within 
this paper.) 

For an open set O C X(0 n dX ^ 0), we call .^-geodesic segments 7 which are 
contained in O with endpoints at dX 0-local £ -geodesies; we denote the set of these 
by Mo- Thus Mo is an open subset of the smooth manifold M of all f -geodesies. 
We then define the local £ -geodesic transform (or the local X-ray transform for a 
general flow) of a function / defined on X as the collection (//) (7) of integrals of 
/ along 7 £ Mo, i-e. as the restriction of the X-ray transform to Mo- We ask the 
following question: 

Can we determine f\o, the restriction of f on the open set O, by knowing 
{If){l),foralljeMo? 
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In order to state our main theorem in conerete terms, we introduce some notation 
below. Let p e C°°(X) be a defining function of dX, so tliat p > Q in X and p < 
on X\X^ vanishes non-degenerately at dX. Our main theorem is an invertibility 
result for the local f -geodesic transform on neighborhoods of p G dX in X of the 
form {x > — c}, c > 0, where a; is a function with x{p) = 0, dx{p) = —dp{p), see 
Section 2 for the definition of x. The statement of our theorem is similar to that of 
[13] on the local geodesic ray transform, we now are able to show the same results 
for the local f -geodesic transform. 

Theorem 1.1. Assume X is strictly £ -convex at p (z dX , there exists a function 
X e C°°{X) vanishing at p with dx{p) = —dp{p) such that for c > sufficiently 
small, and with Op = {x > — c} n X , the local £-geodesic transform is infective on 
H''{Op),s>{). 

Further, let H'^[Aio ) denote the restriction of elements of H'^{Ai) to Mo , 
and for F > let 

HUOp) = e^/^'+^^H%Op) - {/ e HUOp) : e-^/(^+^)/ £ H^Op)}. 

Then for s > there exists C > such that for all f G Hp{Op), 

It is worthy of remark that for this result we only need to assume convexity near 
the point p. This local result is new even in the case that the metric is conformal 
to the Euclidean metric. While this large weight e^/^^"'"^-' means that the control 
over / in terms of If is weak at a; = — c. Here F > can be taken small, but 
non-vanishing. Further, x, whose existence is guaranteed by the theorem, is such 
that a; = — c is concave from the side of Op. 

As an immediate consequence and application of our main theorem, we consider 
a domain X with compact closure equipped with a boundary defining function 
p : X ^ [0,oo) whose level sets p~^{t),0 < i < T for some T > 0, are strictly 
f -convex (viewed from p~^{{t,oo))), and dp is non-zero on these level sets. We 
obtain the following corollary on the global injectivity of the f -geodesic transform, 
which is an analog of the geodesic case in jTS] , the proof is exactly the same. 

Corollary 1.2. For X and p as above, if the set p~^{[T, oo)) has measure, the 
global £- geodesic transform is injective on L'^{X), while if p~^{[T, oo)) has empty 
interior, the global £-geodesic transform is injective on H''{X),s > n/2. 

We also point out that one can derive similar local and global reconstruction 
methods and global stability estimates in the same way as in [13) . 

The (global) X-ray transform has been widely studied, among which the sim- 
pliest case is the geodesic ray transform. The standard X-ray transform, where 
one integrates a function along straight lines, corresponds to the case of Euclidean 
metric and is the basis of medical imaging techniques such as CT and PET. The 
case of integration along more general geodesies arises in geophysical imaging in 
determining the inner structure of the Earth since the speed of elastic waves gener- 
ally increases with depth, thus curving the rays back to the Earth surface. It also 
arises in ultrasound imaging, where the Riemannian metric models the anisotropic 
index of refraction. Uniqueness and stability was shown by Mukhometov [6] on 
simple surfaces, and also for more general families of curves in two dimension. In 
[1], injectivity result was proved for a general family of curves on Finsler surfaces. 
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The case of geodesies was generalized also for simple manifolds to higher dimensions 
in [2], [7], [8]. In dimension n > 3, the paper [3] proves injectivity and stability for 
the X-ray transform integrating over quite a general class of analytic curves with 
analytic weights, assuming an additional microfocal condition that includes the case 
of real-analytic metrics for a class of non-simple manifolds, while [13^ provides a 
completely new approach to the uniqueness for the global geodesic ray transform. 

We remark that Corollary 1.2 provides a new approach to the uniqueness of the 
global problem for the X-ray transform for a general flow. The only method up to 
now, except in the real-analytic category [5] , |10| , has been the use of energy type 
equalities that introduced by Mukhometov [6] and developed by several authors 
which are now called "Pestov identities" . The global geometric condition imposed 
here is a natural analog of the condition -£:{r/c{r)) > proposed by Herglotz [i] 
and Wiechert and Zoeppritz [14 for an isotropic radial sound speed c{r). It allows 
in principle for conjugate points of the flow. In |11| one can find a microfocal study 
of the geodesic X-ray transform with fold caustics. A similar condition of foliating 
by convex hypersurfaces was used in |.l2i to satisfy the pseudoconvexity condition 
needed for Carleman estimats. 

For the local X-ray transform, the only result so far is for the local geodesic ray 
transform proved by Uhlmann and Vasy JL3;. They defined an operator A which 
is essentially a 'microfocal normal operator' for the geodesic ray transform, such 
that Af only depends on the integral of / on the elements oi A4o (the set of Op- 
local geodesies). The operator A is an elliptic pseudodifferential operator only in 
X > — c. However, it turns out that the exponential conjugate Ap of A is a scattering 
pseudodifferential operator in Melrose's scattering calculus fSl on a; > — c, see also 
[T51 Section 2] . They showed that Ap is a Fredholm operator and it is invertible for 
c near 0, which induces both uniqueness and stability estimates for local geodesic 
transform. The key ingredient from the geodesic nature of the curves in their paper 
is that for z € X near point p, and geodesies 7 with 7(0) — z, ^(a;o7)|j^o = 0, then 
Jp-(i o 7)|t=o induces a positive definite quadratic form a near p. This quadratic 
form a plays a crucial role in showing the invertibility of the principal symbol of 
Ap at the boundary x = — c. 

Our main theorem is proved along the lines of [13j . We show that one can prove 
the invertibility of Ap for general flows without the aid of the existence of some 
positive definite quadratic form, by using the polar coordinates to calculate the 
boundary principal symbol of Ap at a; = — c. 

Acknowledgements. The author thanks his advisor Prof. Gunther Uhlmann for 
suggesting this problem and helpful discussions and for reading the previous version 
of this paper. 



2. Proof of the main theorem 

Suppose first that X is a domain in {X,g), p e dX, and dX is strictly £- 
convex at p (hence near p; the convexity assumption will guarantee that in a small 
neighborhood of p, there is no conjugate points with respect to the flow ). That is, 
with p a boundary defining function of X, we have that for verctors v G TpX\o, 

{Hgp){v) = ^ {mp){v) + {E{p,v),Wp{p)) < 0, 
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where (Hgp){v) = (Vp, u), H^p is the Hessian of p. In particular, by the compact- 
ness of the unit sphere and the homogeneity of H^p and E, there is a neighborhood 
J7o of p in X and (5 > 0, Co > such that for vectors v E Tu„X\o, 

\{H,p){v)\ < S\v\, => {HIp){v) + (i?(-,«), Vp) < -CMl 

We want to define a function x near p such that x{p) — 0, the region x > — c, 
p > 0, is compact for c > small, and the level sets of x are concave from the side 
of this region (i.e. the super-level sets of x). By shrinking Uo if needed, we may 
assume that it is a coordinate neighborhood of p. Concretely we let, for e > to 
be decided, and with | • | the Euclidean norm, 

then X > — c gives p + e\z — pp < c and thus p < c; further, with p > this gives 
|z— pp < c/e. Thus for c/e sufficiently small, the region x > — c, p > 0, is compactly 
contained in Ua. Further, for v € TjjgX, Hgx{v) — ~Hgp{v) — eHg\ ■ ~p\'^{v), so 
Hgx{v) = implies \Hgp{v)\ < C"e|w|g, so with (5 > as above there is e' > such 
that for e e (0, e'), Hgx{v) = in t/o implies \Hgp{v)\ < S\v\g, and then for e < e', 

Hp{v) + {E{;v),VS:) 
= -(H^gpiv) + {Ei;v)yp)) e{Hl\ ■ -p\^{v) + {E{;v),V\ ■ -pp)) 

>{C,^C"e)\v\l 
Thus, there is eo > such that for e e (0,eo), H'^x{v) + {E{p,v),Vx{p)) > 
(Co/2)|w|g at TpX when HgX vanishes. Thus taking cq > sufficiently small (cor- 
responding to eo), we have constructed a function x defined on a neighborhood Uq 
of p with concave level sets (from the side of the super- level sets) and such that for 
< c < Co, 

Oc = {x > -c} n {p > 0} 

has compact closure in UoC) X (Here Oc is exactly the Op in Theorem 1.1). 

From now on, we work with Xc = x + c, which is the boundary defining function 
of Xc > 0; we supress the c dependence and simply write x in place of Xc- For most 
of the following discussion we completely ignore the actual boundary p = 0; this will 
only play a role at the end since ellipticity properties only hold in Uq and we need 
/ to be supported in p > 0, ensuring localization, in order to obtain injectivity 
and stability estimates. Thus, completing x to a coordinate system (x, y) on a 
neighborhood Ui C Uq of p, for each point (x, y) we can parameterize 5-geodesics 
through this point by the unit sphere; the relevant ones for us are 'almost tangent' 
to level sets of x, i.e. we are interested in ones with tangent vector k{Xdx + i^dy), 
k{x,y, X,uj) > (to say have unit length), uj E §""^ and A relatively small. 

Now, the iS-geodesic corresponding to (zq, i^q) G SX, 7 = 7zo,i'o' i^ the projection 
of the flow (j)t emanating from (zq, 1^0) (i.e. the intergral curve of Hg through this 
point). If / is a function on the base space X then ^(/ o 7)(0) = Hgf{j'{0)), 

3f(/°7)(0) = HlfiYm + (S(7(0),7'(0)),V/(7(0))) (Here 7 is parameterized 
by arclength). 

We use a fibration by level sets of the function x with non- vanishing differential. 
Letting y be a vector field orthogonal with respect to g to these level sets with Vx = 
1, and using {x = 0} as the initial hypersurface, the flow of V (locally) identifies 
a neighborhood of {x ~ 0} with {—e,e)x x {x = 0}, with the first coordinate 
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being exactly the function x (since time t flow by V changes the value of x to 
t). In particular, choosing coordinates yj on x=0, we obtain coordinates on this 
neighborhood such that dy^ and dx are orthogonal, i.e. the metric is of the form 
f{x, y)dx^+h{x, y, dy) with / > 0. By our convexity assumption, at J^(a;o7)(0) = 0, 
^(xo7)(0) = HlxHm + (i?(7(0),7'(0)),Va:(7(0))) > where 7 = 7.,y,feA.fe^ 
with 7(0) = {x, y), 7'(0) = k{Xdx + tody). 

Now, for given £^-geodesics 7 = ^x,y,kx,kuj (with w e §"^^ and A relatively smaU) 
that are parameterized by arclength, we change the parameters to make the tangent 
vectors 7'(0) = Xd^ + ojdy. So generally the curves 7 are no longer parameterized 
by arclength, we denote the original arclength parameter by s, the new parameter 
by t, and ^ > 0,7^(0) = 7t(0) (i.e. s = 4=> t = 0) . In particular, for 
7 = lx.y.x,u with 7(0) = {x,y), 7'(0) = Xdx+L^dy, if f^{xo-,){0) = (i.e. A = 0, 
thus -±(xo^){0) = too), we have |7'(0)|g = \ujdy\g = 1 (i.e. f |,=o = !)■ By 
convexity assumption, we have 

a(x,y,0,^,0):=i^(xo7)(0) = i^(xo7)(0)>0 

(Roughly speaking, the local change of parameter preserves convexity /positivity) . 
Thus, from now on we use (x, y, A, cu), instead of {x, y, kX, kcu), to parameterize the 
curves 7 for sufficiently small X. 

Similar to |13j . we will work in the following general setting. We consider inte- 
grals along a family of curves ^x,y,x,uj in M", (x, y. A, cj) e R x M"^^ x M x §"^^, de- 
pending smoothly on the parameters. Here K"~^ could be replaced by an arbitrary 
manifold and below we make x small, so we are working in a tubular neighborhood 
of a codimension one submanifold oi X. However, since the changes in the manifold 
setting are essentially just notational (we are working on a single coordinate chart), 
for the sake of clarity we work with M". Further, below we work with neighborhood 
of a compact subset {0} x A' C M^ x MJJ^^^; 7x,y,A,w(0 would only need to be defined 
for (x, y) in a fixed neighborhood U of {0} x K and for |A|, \t\ < Sq, 6q > a. fixed 
constant. 

The basic feature we need is that for a: > and for A sufficiently small, depending 
on X, the curves stay in [0,oo) x M"^^. Thus for a; = 0, only A = is allowed. We 
assume that 

lx,y,x,uj{0) = ix,y), lx,y,\AO) = (A,w), 
7x.,y.\,uji't) = 2(a(a;, y, A, uj, t),l3{x, y, A, cu, t)), 
with a, /3 smooth and 

a(0,y,0,w,0) > 2C > 0. 
This implies that li K C M"^^ is compact, then for a sufficiently small neighborhood 
U of {0} X K in M" (with compact closure), and |A|, \t\ < Sq, where Sq > small, 
we have 

a{x,y,X,uj,t) > C > 0. 
One may assume that x < Sq on U. Thus for j{t) — (x'{t), y'{t)), 

r-l 



X = x + Xt + t^ (1 - C7)a{x, y, A, uj,a) da > x + Xt + Cr /2, 
so if |A|, \t\ < So, {x,y) e U, then 
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Thus, for |A| < \j2Cx (and |A| < (5o), x' > 0, i.e. the curves (f-geodesics) remain 
in the half-space x' > at least for \t\ < Sq. Further, if we fix xq > 0, then x' > xq 
provided |i + ^| > \/2xo/C and \t\ < Sq, thus when |A| < CqXq and |A| < Sq, then 
x' > Xq provided \t\ > ^^^ + y^2xo/C, \t\ < Sq- Assuming x < xq and taking 
Xq sufficiently small so that ^xq + -\/2xo/C < Sq, we thus deduce that the curve 
segments jx,y,\.Lj\{^So,So) ^^^ outside the region x' < xq for t outside a fixed compact 
subinterval of {—Sq, Sq). From now on, by 7 we mean the restriction jx,y,\,uj\{-So.So)> 
and we assume that the functions we integrate along 7 are supported in x' < a;o/2, 
so all integrals are on a fixed compact subinterval. 

Similar to the facts in [3.: and |13| by the convexity assumption, the maps 

T+:SXx [0, <x)^[Xx X; diag], T+{z,v,t) = 7,.,(i) = (z, \z' - z\, ^7^) 

\z'-z\ 

and 

r_ : 5Xx (-^,0] ^ [XxX-diag], r_(z, w,i) = 7,,,(t) = {z,-\z' - z\,^^^—^^ 

\z z\ 

arc two diffeomorphisms near SX x {0}. Since we actually work with (locally, in 
the region of interest) the set of tangent vectors of the form Xd^ + a;9y, where 
w G §"~^ and A relatively small, which can be regarded as a small perturbation of 
a portion of the sphere bundle SX. Thus we reduce 5q if necessary so that r+ is a 
diffeomorphism on Ux,y x {—5o,5{))\ x SJJ^^ x [0,(5o)t, and analogously for r_. We 
assume this from now on. 

Our inverse problem is now that assuming {If){x,y,X,LL>) = J^f{jx,y,\,Ljit))dt 
is known, we would like to recover / from it. Recall our convention from above, the 
integral is really over (— (5o, Sq), and f{x', y') vanishes for x' > a;o/2. One difference 
comparing with the geodesic flow in [T? is that generally the flow defined by E is 
not time reversible, i.e. 

"/x,y,-\,~-ui — t) 7^ Jx,y,\,uj{t). 

Thus we will have two curves with a given tangent line at each (x, y), so having the 
integral of functions along both. 

The idea is similar to [13 , namely to consider for a; > 

(2) {Af){x,y)^ f f {If)ix,y,X,Lo)x{x,X)dXduj, 



where x is supported in |A| < \/2Cx. We define 

X(a;, A) =a;"^x(A/a;), 

where x is compactly supported (for sufficiently small x). We can allow x to 
depend smoothly on y and uj; over compact sets such a behavior is uniform. As 
what mentioned in 13 , for s > 0, 

A = Lol : H'{[0,oo) x M"-i) ^ a;"""^iJ^([0,oo) x R"-^) 

is bounded, where L : H'{[0,oo) x R"-i x M x §"-2) ^ x-''-^H''{[0,oo) x R"^i) 
is also bounded and 

{Lu){x,y) — / / u{x,y, X,Ll!)x{x, X)dXduj. 
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If we show A is invertible as a map between above spaces of functions supported 
near x = 0, we obtain an estimate for / in terms of //. 

Similar to [131 Proposition 3.3], next we show that A is an eUiptic pseudodif- 
ferential operator in x > for proper choice of x; while the conjugates of A by 
exponential weights are scattering pseudodifferential operators on a; > 0. 

Proposition 2.1. Suppose x e C^(R), x > with x > near 0, then A E 
^^^(x > 0) is elliptic, and the operator Ap = x^^e^^^'^Ae^^'^ is in '^^^''^{x > 0) 
for F >0. 

Proof. First we work in a; > 0, ignoring the limit x — ^ 0. The diffeomorphism 
property of T± allows us to rewrite, with \dv\ denoting a smooth measure on the 
transversal such as |(iA||dw|, 

Af{z)^ E f f{z'mr,\z,z'mr,'r{\dv\dt) 

in terms of z, z' as 

(3) / fiz')\z' - z\-^'{biz, \z' - z|, ^^) + b{z, -\z' zl -^^)}dz\ 
J \z - z\ \z' ~ z\ 

b{z, 0, w) — x(z, w)a{z, w), 

where ct > is bounded below (it is actually the Jacobian factor). Thus A is a 
pseudodifferential operator with principal symbol given by the Fourier transform 
of 

z' — z , , _ , , z' — z 



V - zr"+l{(x^)(z, -, r) + (X^)(Z, -, :)} 



in Z ~ z' — z. By similar argument as in the proof of 1131 Proposition 3.3], we can 
get that the principal symbol of A at (z, Q is of the form 

(4) 
CnlCr' / (xa)(z,Z^) + (xa)(z,-Z^)dZ^-2c„|Cr^ / (xa)(z, Z^)dZ^, 

where Z^ is the parameter for the sphere §"~^ which is orthogonal to C,. In par- 
ticular, under our assumption of Xi ^ is an elliptic pseudodifferential operator of 
order —1, in accordance with the results of Stefanov and Uhlmann [9 , provided 
n > 2. 

We now turn to the scattering behavior, i.e. as at least one of x, x' — >■ 0. We 
apply the scattering coordinates (x,y,X, F) introduced in [TJ", Section 2], where 

^^ x' -X ^ ^ y' -y 
x"^ ' X 

With K denoting the Schwartz kernel of A, Ap has Schwartz kernel 

(5) Kp{x, y, X, Y) = x-ie-^^/(i+-^)i^(x, y, X, Y), 

here K has polynomial bounds in terms of X, Y. Our main claim is that Kp and 
its derivatives have exponential decay for F > 0, and Kp is smooth for {X, Y) 
finite, non-zero, conormal to {X,Y) = 0. This will imply that Ap is a scattering 
pseudodifferential operator. 
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Similar to ^^, we can also use {x, y, \y' — y\, \^,_^\ , ttt^) as the local coordinates 
on T+{suppx X [0,(5o)); and analogously for T-{suppx x (— (5o,0]) the coordinates 
are {x, y, —\y' — y\, — i^/I^i , — i^/I^i )■ Now, in terms of the scattering coordinates, 

, , x' — X xX u' — y 
\y' - y\ = x\Yl = — , f^ = r, 

\v -y\ \Y\ \y -y\ 

Thus similar to the estimates in [ISJ by Taylor expansion in \y' — y\, 

xX " 

+ \Y\K{x,y,x\Yl—-,Y), 



A(r+^) _ X ^ „^,x,„ „i^^i xX 



(6) 
Similarly, 



|y| ■ ' ' ^ '^' ' "l^^l 

^ = -^-\Y\Mx,y^-x\Yl-^^-Y). 



xX - 
u;{T^') =Y + x\Y\n{x,y,x\Y\,—,Y), 

tX 

^{r-J) = -Y - x\Y\n{x,y,~x\Yl-—-,~Y) 



and 



(7) 



Thus, 



(8) 



<(r;i) = x\Y\ + x^\Y\^f{x, y, x\Yl — , F), 



xX 

W\ 



tX ~ 

t{r-J) = -x\Y\-x'\Ymx,y,^x\Yl~ — ,~Y). 



dtdXduj{r^^) = J{x,y,±\Y\,±— ±Y)x^\Y\-^dXd\Y\dY 
= J{x,y,±\Y\,±^ ±Y)x'\Y\-"+'dXdY 



where the density factor J is smooth and positive, and J\x=o — 1- Note that on 
the blow-up of the scattering diagonal, {X = 0, F = 0}, in the region |y| > e\X\, 
thus on the support of x in view of ([6]) , 

x,2/,±^,±|r|,±y 

are valid coordinates, with |F| being the defining function of the front face of this 
blow up. Taking into account the x~^ in the definition of x, we deduce that Kp is 
given by 

X tX X 

+x(-^-|>^|A(2:,y,-x|r|,- — ,-f))J(x,y,-|y|,- — ,-y)}, 

so in particular it is conormal to the front face on the blow-up of the scattering 
diagonal, of the form p~"+^6, where b is smooth up to the front face, and without 
the first exponential factor it, together with its derivatives has polynomial growth 
as {X,Y) ^- oo. We decompose Kp into two pieces supported in |(X, y)| < 2 and 
|(X, F)| > 1 by a partition of unity. For the first term, supported in |(X, y)| < 2, 
similar calculations as in ([3]) in Fourier transforming this in {X, Y) show that this 
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term of Kp is indeed the Schwartz kernel of an element of '^^^■^ , with standard 
principal symbol being given by the analogue of Q. So we will show that the 
second term is Schwartz in (X, Y) due to the exponential decay of the first factor 
of ([9]) on the support of x- 
We use ([7]) to express ^ by 

(10) x' ^ X + \t + a{x,y,\uj)f + 0[t^), y' = y + ut + 0{t^), 

where the terms 0{t^) and 0{t^) have coefficients which are smooth in (x, y, \,lo). 

Thus, 

x' — X \t at^ t^ r^, > \ 

X = J— = — + ^- + — JL(a;,2/, A,w,i), 

with T a smooth function of its arguments, so by ([7]) 

X - ^i£ii|r|(i + x\Y\f) + a(r;i)|y|2(i + x\Y\ff + x\y\^t{t^^), 



also 



X = ^i— i|y|(-l - x\Y\f) + a(r:i)|y|2(-l - x\Y\ff - x\Y\'^T{TZ^) 



Thus 



(11) 



A(r-) ^-(r;^)l^^^(,) 



\Y\ 

A(ri^) _ -x + a(r:^)|yp 



+ 0'{x), 



where 0{x) and 0'(a;) have smooth coefficients in terms of a;,?/,x|y|,a;X/|y|, F. 
Thus for ^ e {-c,c), -c\Y\ < X - a{T^^)\Y\'^ < c\Y\, which shows (by the 
positivity of a) that X — > +cx) on suppx if |y| — > oo, and indeed for |y| large 
enough, X > C\Y\'^ for some C > 0. 

Now for all A'' the exponential factor in ([9]) is < C\{X, Y)\^'^ for suitable C on 
the support of x, so combined with the polynomial estimates for the derivatives of 
the other factors, it follows that Kp is smooth in (a;, y) with values on functions 
Schwartz in {X, Y) for {X, Y) ^ 0, and conormal to {X, Y) = 0, which is exactly the 
characterization of the Schwartz kernel of a scattering pseudodifferentail operator. 
This finishes the proof of the Proposition. 

D 

Here the additional information is the behavior of Ap at x = 0, but given that 
Ap is an element of ^jj^'^, the same information can be obtained from computing 
the boundary principal symbol. Indeed this is immediate from Q and (jlip which 
show that at a; = (the scattering front face) the Schwartz kernel K{y,X,Y) of 
Ap is 

^^FX^Y\--+U^^^ X-a{0,y,0,Y)\Y\^ ^ ^ ^^ -X + a{0,y,0,-Y)\Y\^ ^^^_ 

So the desired invertibility of Ap amounts to the Fourier transform of the kernel, 
K{y, •, •) being bounded below in absolute value by c((f , 77))^^, c > (here {£_, 77) are 
the Fourier dual variables of {X,Y)). Thus our job now is to compute the Fourier 
transform of K{y,-,-). Note that we allow that x also depends on y and w. We 
have thus shown 
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Lemma 2.2. The boundary principal symbol of x^^e^^'^Ae^'^ is the {X^Y)- 
Fourier transform of 

"" +.( -'""°'°,^°-'^''"^l' .,.-f-)). 

In order to find a suitable x to make Ap invertible, we follow the same strategy 
of [13], namely we do calculation for x(s) = e~* ^^^^ "■* with F > { here we need 
the positivity of a), so x(-) = c\/F~^ae~^ "'I /^ for appropriate c > 0. Here x 
does not have compact support, and an approximation argument will be necessary. 
First, the Fourier transform of -ftT in X is 

where a+ — a{0,y,0,Y), a- = q;(0, y, 0, —y). Substituting the particular x into 
(|13p yields a non-zero multiple of 

(14) |r|2-"{(F-ia+)^e-^"(«'+^')"+l^l'/2 + (F-V_)^e-^"(«'+^')"-l^l'/2}. 

Now we estimate the Fourier transform of (|14p in y upto some non-zero multiple. 
Not like the case of the geodesic flow in [13], generally our a may contain terms 
other than a quadratic form in Y, thus we use polar coordiantes to compute the 
Fourier transform in Y. We denote '^^ — - by 6, then 

^y(^xAl(y,C,^) - /"e-^''-^|r|2-"(ay'e-''"+l^l' +aL/'e-^"-l^l')dy 



+ 00 



2 Js,i-2 

^5-1/2 f p-\rrY\''/'ib»{y,Y)^Y^ 

JS"-2 

here ~ means equal up to some multiple. 

To estimate the Fourier transform of (IT^ . we need to study the joint {£,,ri)- 
behavior, i.e. when {{^,r])) is going to infinity, where we need lower bounds, we 
denote {^^ j^ p'^y/'^ by (^)^ then J^x.YK{y,S,,'q) is a constant multiple of 

(0-^ / ^-\Tb-y\'/^-<v^y)dY 

with c = c{F) > 0. The only property of a{y, Y) we need is locally uniformly 
positivity, i.e. < ci < a < C2 for some positive constants ci,C2 that depend on y 
and are locally uniform, which is true under our assumption. 
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When M < C, then (S,)^^ is equivalent to ((^, »7))^^ in this region in terms of 
decay rates, 

^-\^.Y\^„c.iy.Y)^Y> I e~''\Tb\'dY> j e-~^dY = C. 

Thus J-x,yK(y,e,r/) > C{ir^ ~ C{{(„7^))~\ 

When jA is bounded from below, in which case ((^, r/))~^ is equivalent to |»y|~^, 
we write F = (F H , Y^) according to the orthogonal decomposition of Y relative to 
1^, where Y^ =Y ■ ^, and dY is of the form a{Y\\)dY\\de with 9 = ^, a(0) = 1 
then 



e-='(^"lt|)^a(yii)d0dyii 



(15) 






(0 


bfr'"*' 


}a(fll)dfl 


'I 


de. 

-3 








Since 




-c'(yII 


I>7| -12 

to) ^,5c 


1 in distributions as 


, ^^^ -> oo, 


1151 


is equal to 


hi 


/s^ 


-3 de = 


2^lf} 


(C: 


> 0) modulo terms decaying faster 


as M^c^. Since /JM, 


^ — C 


:'(yii 


hi -,2 - II - II 

to) }a(yii)dyii 


is smooth 


inM 


, there is 


iV > 0, such that 















{Me-'(^"w'^Hyii)dyii / de>c 

for jA > N . (Notice that the integral on S"^'^ uses very strongly the assumption 
n > 3; when n — 3, d0 is the point measure) Thus Tx.YK{y,^,ri) > C't^tH = 

Thus we deduce that J^x,YK{y,£_,ri) > c{{^,i])) ^ for some c > 0, i.e the ellip- 
ticity claim for the case that x is a Gaussian. By approximation argument as in 
[13] , one obtain some x € C*^ (R) such that the Fourier transform of K with this x 
still has lower bounds c{{^,ri))^^,c > 0, as desired. We have thus proved; 

Proposition 2.3. For F > there exists x G C^{M.), x > 0, x(0) = 1, such 
that for the corresponding operator x^^e^^'^Ae^'^ the boundary principal symbol 
is elliptic. 



Ap = x-'e-^/^Ae^/- e ^r/'" 



Hence, we have 

Af 

is elliptic both in the sense of the standard principal symbol (in the set of in- 
terest Oc), and the scattering principal symbol, which is at a; = 0. In particu- 
lar the results of [13l Section 2] on elliptic scattering pseudodifferential operators 
are applicable. Thus, for c smaU, if we denote Ale — {xc > 0}, H^'^{Mc)k = 
{/ S Hgi!'{Mc) ■ suppf (- K {K is a compact subset of Oc)} (Roughly speak- 
ing, Hli!'{Mc) = x'^^H''{Mc), see [13l Section 2] for the detailed definition and 
properties), 
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satisfies the estimates 

ll/lle-f/-//|i'-(M<:)A- < C'P/lleJ=-/-ff| + i-'-+i(Me)- 

By similar arguement as in the last part of [T3], we thus deduce for s > 0,6 > 

|!/||e(F+5)/xfl-.-l(jif^)^^ < C\\Af\\eP/^H-{M,) < C!' W^ f\\ H- (Mm^) 

when / e H''{Mc)k- With F + (5 replaced by F (since both F > and 5 > are 
arbitrary), this completes the proof of the main theorem. Thus we obtain the local 
injectivity and stability estimates for the X-ray transform for general flows. 
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